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We show that the 3 + 1 vacuum Einstein field equations 
in Ashtekar's variables constitutes a first order symmetric hy- 
perbolic system for arbitrary but fixed lapse and shift fields, 
by suitable adding to the system terms proportional to the 
constraint equations. 



I. INTRODUCTION 

In the mid-eighties a new representation for Hamilto- 
nian General Relativity was introduced by Ashtekar [Q. 
The Ashtekar representation maybe understood as result- 
ing from a complex canonical transformation that goes 
from the ADM gravitational phase space variables to a 
new set of variables. One of the major benefits of this 
representation is that the constraints take a simpler form 
in terms of these variables than in the ADM formulation. 
Another important aspect of this set of variables, surpris- 
ingly not yet exploited and not shared by the standard 
ADM representation, is that the evolution equation are a 
first order system of differential equations. Thus it is nat- 
ural to consider the well-posedness of the classical initial 
value problem in this context. 

Well poscdncss of Einstein equations in itself is a long 
ago closed issue. Since the initial data for the equations in 
Ashtekar's variables is in one to one, and continuous cor- 
respondence with the ADM variables, the result follows 
at once. Rather, the issue we pose here is whether Ein- 
stein's equations in Ashtekar's variables conform a good 
system to evolve arbitrary initial data using approxima- 
tions schemes, be numerical or analytical. Bad choices 
of evolution gauges usually lead into spurious coordinate 
singularities which effectively destroy the smoothness of 
the map between variables. We are also concerned with 
consistency, convergence and stability properties of ap- 
proximation methods which use these variables. These 
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issues have not been fully solved in the realm of the nu- 
merical analysis of non linear partial differential evolu- 
tion equations, but all powerful results so far obtained 
make heavy use of the symmetric hyperbolic character of 
the equations. This is the main reason for checking this 
property here. Updated references to this problem are in 
§, |, §, [|, and [|. 

As discussed in the conclusions (section |y|), one im- 
portant property of the system we consider is that the 
lapse-shift pair can be given before hand on spacetime or 
can be obtained along evolution, either imposing evolu- 
tion equations and so effectively enlarging the symmetric 
hyperbolic system, or imposing elliptic equations on each 
hypersurface t = const. 

The plan of the paper is as follows. Section JO] is 
devoted to mathematical preliminaries where we recall 
some aspects of Ashtekar's Hamiltonian formulation. In 
section 



III we present the main result. 



II. THE FIELD EQUATIONS IN ASHTEKAR'S 
VARIABLES 



In this section we shall present the necessary defini- 
tions and a brief overview of the Hamiltonian formalism 
in spinorial variables in order to get Einstein field equa- 
tions as a system of evolution and constraint equations 
(we follow §, |, and §). 

Let M be a four dimensional manifold with a smooth 
metric g a b of signature ( — h H — h). Furthermore, assume 
that M admits a time function t, a smooth scalar field on 
M, whose gradient is everywhere timelike and denote the 
time direction by t a , which is a smooth vector field on 
M with affine parameter t, i.e. t a V a t = 1. Then the one 
parameter family of surfaces E t , defined by t = const., 
are three dimensional spacelike surfaces. Finally we shall 
denote by n a the future directed unit normal vector field 
to E t , and the intrinsic and positive definite metric on 
Et by q a b- By projecting t a into and orthogonally to E t , 
we obtain the lapse and the shift fields N and N a , thus 

We shall also assume that our spacetime admits an 
spinorial structure (see |ioj , and references therein ); i.e 
we shall assume that to each point p E M we can assign a 
complex two dimensional vector space W equipped with 
a nondegenerated symplectic form (skew 2-form) . The el- 
ements of W will be denoted by £ A and the corresponding 
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symplectic form by cab (its inverse will be denoted by 
e AB , and is such that eac£ AB = 5c B )- The complex con- 
jugate vector space associated with W will be denoted 

_A' 

by W, its elements by £ and the symplectic form by 
£A' B'fjJ Since the group that preserves the structure on 
W is SL(2,C), £ A will be called a SL(2,C) spinor at 
p G M. 

To relate tensors to spinors on M, one must fix a 
metric preserving isomorphism between the four dimen- 
sional real vector space V = {£ AA € W <g> W : ^ Aj4 ' = 
— i3aa>bb< = tAB^A' b' } and the tangent space 
T p (M) with metric g a f,. This isomorphism is called an 
SL(2, C) soldering form and is denoted by <J a A A>, Then 

ah a b AB A' B' 

g =cr AA'cr BB>£ e 

The covariant derivative operator on tensors has a unique 
extension to spinors if we impose the condition 



A A' 



= 0. 



In order to relate spinors on these surfaces to their in- 
trinsic geometry, we shall need an additional structure 
on W which reduces the structure group SL(2,C) to 
SU(2). This additional structure on W is a positive 
definite hcrmitian inner product which is denoted by 
G : W x W -> C, or 

Gaa'-.WxW^C, (r) A ,f Hf'W- 

In this notation hermiticity appears as Ga'A — Gaa 1 - 
We will assume that (either cab or) Gaa< is so normal- 
ized that 

e A B Ga'aGb'b = £ab, or Ga>aG a b — 5a b , 

In the formalism 3+1 this additional structure is obtained 
by choosing 



G 



AA> 



then as a a AA = —o~ a AA and n a is timelike and future 
directed, we have the positive definite hermitian inner 
product that we were looking for. Via this inner product, 
we identify unprimed SL(2, C) spinors on M with SU(2) 
spinors on £ t . The relationship between the intrinsic 
geometry with SU (2) spinors is given by: 

_ A . „ ft AA' n 
&a B ■— <la &b ^A'B 

that is the SU(2) soldering form. Then 



V2i 



nAA> 



V2 



AA' 



q a b 



-a a A <JbB 



-tr(a a a b ). 



Here a matrix notation is employed (and shall be used in 
the following) for unprimed spinor indices in which ad- 
jacent summed indices go from upper left to lower right, 
e.g., {a a a b ) A B =a a A C <J b c B . 

There is a natural (canonical) spinorial connection as- 
sociated with the metric such that D a a b = 0. The pull- 
back of the 4-dimcnsional canonical connection to the hy- 
persurface is another connection on the 3-surface called 
Sen Connection, defined by 

V a a A = q a b ^7 b a A , 

the curvature associated to this connection is denoted by 
FabA B and is the pullback of the 4-dimensional curvature 
F a bA B = q a c qb d i RcdA B - If we fix (for simplicity) a 
connection d on the unprimed spinor indices, chosen to 
be real and flat, and replace V by a SL(2, C) Lie algebra- 
valued connection 1-form 4 A a A B - 

^a^A ■= d a ^A + 4 A aA B ^B, 

the curvature tensor becomes 

*R ab :=2d [a 4 A b] + [ 4 A a , i Ab}. 

Recall that 4 R ab A B is defined by 

(V Q V h - V b V a )X A = 4 R abA B \ B , 

thus we can write F ab c D in terms of a SU(2) Lie algebra 
valued one form A a c D 

F ab := 2d [a A b] + [A a ,A b ]. 

Using these variables the action can be defined as: 

S(o; 4 A) = J d 4 x ( 4 o-)a a A A '<J b BA> 4 Ra b AB , 

where ( 4 <r) denotes the square root of minus the deter- 
minant of the four metric, this action is presented in 0] 
and in jjj . Using the identity: 

a A'b 4p AB ,„ / , /o „• „abi p „a Mp I 

a a o BA> Rab = tr I \J i i n a H ab — a a H ab I 

= trf- cr a o- bi R ab + V2 i ArV 
x{C t 4 A b -V b { 4 A-t)) 
— \[2 i N~ 1 N a a b4 R ab ^J , 

where 4 A ■ t — 4 A a t a , we have 
S = J At J d?x( A a) trU/2 iN- 1 o- b (£ t 4 'A b 

-V b { 4 A ■ t)) - G a a h4 R ab - V2 iN~ 1 N a a b4 R ab ^ . 

Finally, since a a ab projects onto the 3-surface, and 4 ct = 
Na, the actionn becomes 
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S = J At J d 3 x tr (V2 i& b (£tA b -VX A -t)) 

-N d a d b F ab - V2 iN a d b F ab ^j . 

with a a AB = (cr) a a AB and N = (cr) -1 AT. 

The action depends on five variables N, iV a , 4 A • t, 

A aA B and o a ab • The first three variables play the role 
of the Lagrange multipliers, only the last two are dynam- 
ical variables. Varying the action with respect to the 
Lagrange multipliers we obtain the constraint equations: 



C(a,A) 
C a (a,A) 
C A B (a,A) 



tr{a a a b F ab ) = 0, 
tr(t b F ab ) = 0, 
V a a a A B = 0. 



(1) 



varying with respect to the dynamical variables, yields 
the evolution equations: 

C t a b = [ 4 A ■ t, d b ] + 2V a (N^ a ^) 

C t A b = V b { 4 A ■ t) + N a F ab + ±=i N[a a , F ba ]. (2) 



III. SYMMETRIC HYPERBOLICITY 

In order to obtain a symmetric hyperbolic system for 
any given lapse-shift pair, we suitably modify the field 
equations outside the constraint submanifold. Indeed, 
the modification only involves the addition of terms pro- 
portional to the constraints. Since we know that the evo- 
lution vector field is tangent to this manifold, the physi- 
cal evolution, that is the dynamics inside the constraint 
manifold, remains unchanged. The modified evolution 
equations are 

£ t a b = [ A A ■ t,a b ] + 2V a {N^ a ^) 

--^W a (N[d a , a b }) + Np, ~a b \ + N b C, 
V 2 v 2 

C t A b = V b ( A A ■ t) + N a F ab + -Li N[a a , F ba ] 



' m b c + — A Wb dc ~PcC d . 



>V2 



(3) 



Since the Lagrange multiplier 4 A a t a is an arbitrary func- 
tion we can choose it as 4 A a t a = A a N a , i.e. we set its 
normal projection to zero. In this way, in the equation 
for the evolution of A a , the principal part correspond- 
ing to the first two terms is replaced by the directional 
derivative of A a . 



To see that the system (||) is symmetric hyperbolic, 
we shall only need to consider its principal symbol. Re- 
call that the principal symbol of a quasilinear evolution 
equation system 

it 1 = B l j a (u)V a u j + M l (u) 

is given by P(u, ik a ) = i B l j a (u)k a . In our case u denotes 
u=(a a ,A b ). 

Lemma III.l: The equation system (^j for any fixed, 
but arbitrary lapse and shift fields is a symmetric hyper- 
bolic system. 
Proof: 

Since the vector u in our case is complex, symmetry of 
real systems corresponds to antihermiticity of the prin- 
cipal symbol with respect to some local hermitian inner 
product, i.e. an hermitian, positive definite bilinear form. 
In our case the natural one: 

(«2,«i) = ((a 2 ,^ 2 ),(a 1 ,A 1 )> 

= tr(ert"^)g a6 +tr(^ 2 t A 1 6 ) q ab , 

is the one that works. Thus we need to prove that 

P12 + P^2 = (u 2 , + (tta, PW) = 0. 

To obtain an expression for the symbol Pu± we neglect all 
the terms which do not have derivatives of the dynamical 
variables a b and A b . We remark the very important fact 
that for both equations, (||) and (||), the principal part 
decouples into two blocks one acting on a b and the other 
acting on A ai so we can write 



P(u,ik a ) = P a (u,ik a ) © P A (u,ik a ), 



(4) 



in our casey = (a, A). The action of P on any solution 
u can be written as 

P a (&, ik a )ai b = L= N[a a , a x b ] + i k a N a ^ b , 

P A (a, ik c ) A 1 a = i N b k b A\ - -4=N [a b , k a A\] 

v2 

+i=N [d b ,k b A\] 

-i=N a a tr([a d ,a b ]k d A\) (5) 
-N e a bc a c ix(a*k b A 1 d ) 
+N i a bc & c tx{a d k d A\). 

Adding the expressions for P a n and PJ12 we obtain 



= -^=mr(a a [a 1 b ,a 2b }) 
+^|Ntr(a> lb ,a 2 b ]) 



(6) 



ik a N a tiia 2b a 1 b ) + ik a N a ti(a2 b a lb ) 



0. 
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here we have used that a a > — —a a . It only remains to 
prove a similar result for Pa 12, using in the expression 
for Pa A 1 a that 



and 



we obtain 
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Pa A 1 a 



A 1 a = -tY^ad^d 6 



\^.a ~b] B 



V / 2 ~ ab -c B 



~bem 



a m k a tr(A b a e ) 



N 



N k b tv(A\a e )a 



~bem ~ 



a m k b tr^^CTe) - e dbe a a k d tr(AVe) 
-~e a bc a c k b tr(AV d ) + e a bc a c k d tr(A 1 6 cr d ). 



very practical, for one would like, for instance, to choose 
the pair along the evolution in order to avoid the occur- 
rence of spurious coordinate singularities. There are at 
least two ways to loosen the assumption made. The first 
is to incorporate the lapse-shift pair into the evolution 
equations by choosing appropriate equations for them in 
such a way as to obtain a bigger symmetric hyperbolic 
system. The second is to impose some elliptic equations 
on the pair in such a way as to be able to estimate the 
L norm of the derivatives of it -for they appear on the 
equations for (a, A)- in terms of the L 2 norms of (a, A). 
This would be the case if one were to use Witten's equa- 
tion to determine the pair. Besides these two ways of 
fixing the lapse-shift pairs along evolution there does not 
seem to be any other obvious candidate. 



Thus we have 
N 



2 

(PA12 + P f Ai2) = -V N " k b tv(A 2 ^a e )tv(A\a e ) 



N 



+ ^ N b k b tr(A 2 K^ e MA la a e ) 



+e bem ti{A 2U ° m )k a tr(AVe) 



— e 



■bem 



tr(A 2 ^ a d m )k b tr{A 1 a a e ) 
-i dbe tr(A 2 ^a a )k d tr(AVe) 
tr(A 2 ^ a a c )k b tr(AV d ) 



— e 



- be 



+e a bc tr(A 2 ^~a c )k d tr(A\d d ) (7) 

+ i bern ^2^-^ t r(,4 lo * m ) 
_ ? 6em t T (A la a m ) 

~~e dbe tr(A 2 \d e )k d tr(A la a a ) 
tv{A^ d d d )k b tv{A la d c ) 



-e, 



be 



+ ~e a bc tr(A 2 \v d )k d tr(A la a c ). 

Rearranging the indices and using the antisymmetry of 
the Levi-Civita tensor density we get that the first term 
cancels out the second, the Ath cancells out the 9th , and 
the 5th and 6th cancel 10th and 11th. 
Finally, using 



Ozlf a ^.^l ~abe~ Ad~c 

IA 1 a ' = — e e dce A a , 



[1] New Hamiltonian formulation of general relativity, A. 

Ashtckar. Physical Review D, 36, 6, 1587, 1987. 
[2] Hyperbolic reductions for Einstein's equations, H. 

Friedrich, Class. Quantum Grav., 13, 1451-1469, 1996. 
[3] Geometrical Hyperbolic systems for General relativity and 

Gauge theories, A. Abraham-A. Anderson-Y. Choqut 

Bruhat-J. York, Class. Quantun Grav., 14, A9-A22, 

1997. 

[4] New Formalism for Numerical Relativity, C. Bona-J. 

Masso-E. Seidel-J. Stela, Physical Review Letters, 75, 

4, 600-603, 1995. 
[5] Horizon Boundary condition for black hole spacetime, 

P. Anninos-G. Daues-J. Masso-E. Seidel-Wai-Mo Suen, 

Physical Review D, 51, 10, 5562-5578, 1995. 
[6] First-Order Symmetric Hyperbolic Einstein's Equations 

with Arbitrary Fixed Gauge, S. Fritelli and O. reula, 

Physic. Rev. Letters, 76, 4667-4670, 1996. 
[7] New Perspectives in Canonical Gravity, A. Ashtekar, 

edited by Bibliopolis, 1988. 
[8] Quantum Theory of Spin-?,/ '2 Fields in Einstein Spaces, 

A. Sen. International Journal of Theoretical Physics, 21, 

1, 1982. 

[9] Covariant Action for Ashtekar's Form of Canonical 
Gravity, T. Jacobson and L. Smolin, Class. Quantum 
Grav., 5, 583-594, 1988. 
[10] Spinors and space-time, R. Penrose and W. Rindler, 
edited by Cambridge University Press, 1, 1984. 



the remaining terms vanish, yielding the antihermiticity 
of the principal symbol. This concludes the proof of the 
Lemma. 4k 



IV. CONCLUDING REMARKS 



The result obtained here assumes, for simplicity, a 
jiven and fixed, lapse-shift pair. This assumption is not 
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